Excavation method by blasting

ABSTRACT

The present invention relates to a blasting method which comprises conducting a delay blast at a particular location; predicting time series data of a waveform of ground vibration or noise at a remote location to be generated by a hypothetical single-hole blast at the particular location using at least one of previous time series data of a waveform of ground vibration or noise generated by said delay blast and actually monitored at the remote location, and the corresponding previous actually applied initiation time series of said delay blast; computing a delay blasting initiation time series for a delay blasting, which provides a waveform of ground vibration or noise satisfying specific conditions, based on the above-predicted time series data of a single-hole blast; and carrying out a subsequent delay blast according to the computed delay blasting initiation time series.

TECHNICAL FIELD

The present invention relates to a blasting method capable of reducingground vibration and noise generated upon blasting.

BACKGROUND ART

Conventionally, delay blasting methods using a delay detonator have beenmost advantageously employed to reduce ground vibration or noiseeffectively upon blasting. As methods for reducing ground vibration ornoise more effectively, Japanese Patent Publication No. 122559/1995,Japanese Patent Application Laid-Open No. 285800/1989 and the like haveproposed blasting methods using a detonator excellent in time accuracywhich is controlled by integrated circuits, wherein dominant frequencyor a waveform generated by a test single-hole blast is preliminarilymonitored at a location where the ground vibration or noise becomesproblematical and initiation intervals for a delay blast are determinedbased on the above-monitored dominant frequency or waveforms.

The waveforms of the ground vibration or noise generated by a blast aregreatly influenced by the type of a target rock. In order to reduceground vibration or noise generated by blasting a target rock mosteffectively according to the above methods, it is necessary to monitordominant frequency or waveform of ground vibration or noise which isgenerated by a test single-hole blast at problematic locations everytime before blasting a target rock.

Therefore, it is difficult to minimize ground vibration or noiseconstantly according to the conventional methods.

DISCLOSURE OF THE INVENTION

For avoiding the above drawback, the present invention provides ablasting method comprising conducting a delay blast at a particularlocation; predicting time series data of a waveform of ground vibrationor noise at a remote location to be generated by a hypotheticalsingle-hole blast at the particular location using at least one ofprevious time series data of a waveform of ground vibration or noisegenerated by said delay blast and actually monitored at the remotelocation, and the corresponding previous actually applied initiationtime series of said delay blast; computing a delay blasting initiationtime series for a delay blasting, which provides a waveform of groundvibration or noise satisfying specific conditions, based on theabove-predicted time series data of a single-hole blast; and carryingout a subsequent delay blast according to the computed delay blastinginitiation time series.

The present invention relates particularly to a blasting methodcomprising conducting a delay blast at the particular location; thencomputing the Fourier Transform of the time series data of a waveform ofground vibration or noise generated by said delay blast and actuallymonitored at the remote location, and the corresponding actually appliedinitiation time series data of said delay blast to obtain correspondingspectrums; predicting spectrums corresponding to time series data of awaveform of ground vibration or noise at a remote location to begenerated by a hypothetical single-hole blast at the particular locationusing the corresponding spectrums obtained in the previous step;performing with the spectrums; computing the Inverse Fourier Transformof the performed spectrum; predicting time series data of a waveform ofground vibration or noise at the remote location to be generated by saidhypothetical single-hole blast at the particular location; computing adelay blasting initiation time series for a delay blasting, whichprovides a waveform of ground vibration or noise satisfying specificconditions, based on the above-predicted time series data of asingle-hole blast; and carrying out a subsequent delay blast accordingto the computed delay blasting initiation time series.

The present invention also relates particularly to a blasting methodcomprising conducting a delay blast at the particular location; thencomputing the cross-correlation sequence of time series data of awaveform of ground vibration or noise generated by said delay blast andactually monitored at a remote location, and the auto-correlationsequence of the corresponding actually applied initiation time seriesdata of said delay blast; predicting time series data of a waveform ofground vibration or noise at a remote location to be generated by ahypothetical single-hole blast at the particular location, which mostcertainly seems to form the time series data of a waveform of groundvibration or noise of said delay blast, by solving Wiener's leastsquares theory according to the Levinson algorithm; computing a delayblasting initiation time series for a delay blasting, which provides awaveform of ground vibration or noise satisfying specific conditions,based on the above-predicted time series data of a single-hole blast;and carrying out a subsequent delay blast according to the computeddelay blasting initiation time series.

It is possible to exemplify various methods for predicting time seriesdata of a waveform of ground vibration or noise at a remote location,which is to be generated by a single-hole blast, using time series dataof a waveform of ground vibration or noise generated by a delay blast ata particular location and the delay blasting initiation time series ofsaid blast. The present invention may employ either a method which onlyuses the ground vibration or noise time series of a current delay blast,i.e,, a latest delay blast, and delay blasting initiation time series ofsaid blast; or a method which uses the time series data of groundvibrations or noises of several previous delay blasts besides thecurrent delay blast and delay blasting initiation time series of saidprevious blasts. In order to provide a clearer idea on the presentinvention, there will be described hereinafter several examples of themethod which employs only the time series data of ground vibration ornoise of a current delay blast and delay blasting initiation time seriesof said blast.

First of all, a successive analytical prediction method is described.

Defining the time series data of ground vibration or noise generated bya current delay blast at a particular location and a delay blastinginitiation time series of the blast as a_(m) and Δ_(i), respectively,the time series data X_(m) of ground vibration or noise generated by asingle-hole blast to be predicted can be successively computed as shownbelow. Both a_(m) and X_(m) indicate an m^(th) data sampled under theconditions of a sampling interval of Δ_(t) and a number of samples of N.Accordingly, m falls within the range of 0≦m≦N−1. Δ_(i) is an integerobtained by dividing i^(th) delay blast initiation time T_(i) withΔ_(t). When the number of periods is defined as L, i falls within therange of 0≦m≦L−1. In this case, Δ₀ indicates 0.

Δ₀ ≦t≦Δ ₁ , X _(t) =α _(t)

Δ₁ ≦t≦Δ ₂ , X _(t) =α _(t) −X _((t−Δ1))

Δ₂ ≦t≦Δ ₃ , X _(t) =α _(t) −X _((t−Δ2)) −X _((t−Δ2))

Δ_(i) ≦t≦Δ _(i+1),  Δ_(L−1) ≦t≦N−1,$X_{t} = {a_{t} - {\sum\limits_{n = 1}^{i}X_{({t - {\Delta \quad n}})}}}$

$X_{t} = {a_{t} - {\sum\limits_{n = 1}^{L - 1}X_{({t - {\Delta \quad n}})}}}$

Next, the Fourier Transform method is described.

Defining the time series data of ground vibration or noise generated ata particular location by a current delay blast as A_((t)), delay blasttime series data of the blast as ζ_((t)), and time series data of groundvibration or noise of a single-hole blast to be predicted as X_((t)),the following relationship is recognized among the three kinds of timeseries data.$A_{(t)} = {{\sum\limits_{s = 1}^{n}{X_{({t - {ts}})} \cdot \zeta_{({ts})}}} = {X_{(t)}*{\zeta_{(t)}{()}}}}$

Namely, the waveform A_((t)) derived from a delay blast is representedby a convolution of the waveforms X_((t)) of a single-hole blast andζ_((t)), wherein t₀=0 and X_((t))=0 when t<0.

Supposing, for example, the amplitude of each period is the same,ζ_((t)) becomes 1 when an initiation timing t is t₀, t₁, . . . andt_(n), and it becomes 0 when t is other than t₀, t₁, . . . or t_(n).

Computing the Fourier Transform of the above equation:

A _((f)) =X _((f))·ζ_((f))(f: Frequency)

Accordingly,

X _((f)) =A _((f))/ζ_((f))

Since A_((f)) and ζ_((f)) are known from A_((t)) and ζ_((t)), X_((f)) isobtained. The next steps comprises computing Inverse Fourier Transformof the thus-obtained X_((f)) in order to transform X_((f)) from afrequency region to a time region and obtaining time series data X_((t))of ground vibration or noise of a target single-hole blast to bepredicted.

Next, the de-convolution method is described.

Defining the time series data of ground vibration or noise generated bya current delay blast at a particular location as A_(t), ideal groundvibration or noise time series data obtained by eliminating errors ofmeasurement and correlating deviation among each single-hole blast asB_(t), delay blast initiation time series data of the blast as ζ_(t)(supposing the amplitude of each period is the same, ζ_(t) becomes 1when an initiation timing t is t₀, t₁, . . . and t_(n) and it becomes 0when t is other than t₀, t₁, . . . or t_(n)), and time series data ofground vibration or noise of a single-hole blast to be predicted asX_(t), the following relationship is recognized among the four kinds oftime series data.${\sum\limits_{s = 0}^{m}{X_{s} \cdot \zeta_{t - s}}} = {{X_{t}*\zeta_{t}} = {B_{t} \approx {A_{t}\left( {\,^{*}{:{Convolution}}} \right)}}}$

If it is possible to compute X_(t) so as to make the error between A_(t)and B_(t) minimum, the computed X_(t) will be the ground vibration ornoise data of a single-hole blast to be intended to obtain.

The ground vibration or noise data of a single-hole blast is obtained inaccordance with the following method according to Wiener's least squarestheory.

First, defining the energy of the error between A_(t) and B_(t) as E,the following equation can be established.$E = {\sum\limits_{t = 0}^{n}\left( {A_{t} - B_{t}} \right)^{2}}$

Further, $B_{t} = {\sum\limits_{s = 0}^{m}{X_{s} \cdot \zeta_{t - s}}}$

Consequently,$E = {\sum\limits_{t = 0}^{n}\left( {A_{t} - {\sum\limits_{s = 0}^{m}{X_{s} \cdot \zeta_{t - s}}}} \right)^{2}}$

The energy of the error becomes minimum when ∂E/∂X_(i)=0. Therefore,$\begin{matrix}{{{\partial E}/{\partial X_{i}}} = \quad {{\partial\left\{ {{\left( {\sum\limits_{t = 0}^{n}A_{t}} \right)2} - {2{\sum\limits_{t = 0}^{n}{A_{t}{\sum\limits_{s = 0}^{m}{X_{s} \cdot \zeta_{t - s}}}}}} + {\sum\limits_{t = 0}^{n}{\left( {\sum\limits_{s = 0}^{m}{X_{s} \cdot \zeta_{t - s}}} \right)2}}} \right\}}/{\partial X_{i}}}} \\{= \quad {{{{- 2}{\sum\limits_{t = 0}^{n}{A_{t}\zeta_{t - i}}}} + {2{\sum\limits_{t = 0}^{n}{\left( {\sum\limits_{s = 0}^{m}{X_{s} \cdot \zeta_{t - s}}} \right)\zeta_{t - i}}}}} = 0}}\end{matrix}$

Accordingly,${\sum\limits_{s = 0}^{m}{X_{s}{\sum\limits_{t = 0}^{n}{\zeta_{t - s}\zeta_{t - i}}}}} = {\sum\limits_{t = 0}^{n}{A_{t}\zeta_{t - i}}}$

wherein${\sum\limits_{t = 0}^{n}{\zeta_{t - s}\zeta_{t - i}}} = {\varphi_{i - s}\left( {\varphi:{{Auto}\text{-}{correlation}\quad {function}\quad {of}\quad \zeta}} \right)}$${\sum\limits_{t = 0}^{n}{A_{t}\zeta_{t - i}}} = {\psi_{i}\left( {\psi:{{Cross}\text{-}{correlation}\quad {function}\quad {of}\quad A\quad {and}\quad \zeta}} \right)}$

Consequently,${\sum\limits_{s = 0}^{m}{X_{s}\varphi_{i - s}}} = \psi_{i}$

The aimed waveform X_(t) formed by a single-hole blast is computed bysolving the above equation according to the Levinson algorithm.

In order, to make more precise predictions according to these methods,it is necessary to make a SN ratio of time series data obtained by acurrent delay blast at a particular location as good as possible using adisplacement averaging, a band pass filtering and the like.

Further, there can be suggested several methods for computing, based onthe above predicted data of a single-hole blast, a delay blastinginitiation time series forming a waveform of ground vibration or noiseof the delay blast which satisfies specific conditions. For example,there is exemplified a method disclosed in Japanese Patent PublicationNo. 122559/1995 wherein initiation time intervals are set based on thedominant frequency so as for a wave to interfere with each other; amethod disclosed in Japanese Patent Application Laid-Open No.285800/1989 wherein waveform of the blast is predicted based on thesuperposition theorem to select an optimum time interval; a methoddisclosed in Japanese Patent Publication No. 14480/1996 wherein M seriesis used; a method disclosed in the Journal of the Japan ExplosivesSociety, NIPPON KAYAKU GAKKAI-SHI, vol. 55, no. 4, 1994 whereinauto-correlation and cross-correlation functions are used; and the like.

The specific conditions mean to minimize evaluated values such asdisplacement amplitude, displacement velocity amplitude, displacementacceleration amplitude, vibration level, vibration acceleration level orthe like in the case of a wave, and to minimize evaluated values such assound pressure amplitude, noise level or the like in the case of anoise. Sometimes, the specific conditions mean to minimize the aboveevaluated values in the specific range of frequency.

Once the delay blasting initiation time series is computed, a blast iseffected according to the computed time series with a detonatorexcellent in time accuracy which is disclosed in, for example, JapanesePatent Application Laid-Open Nos. 261900/1987 and 285800/1989. Theground vibration or noise derived from the blast is monitored at aspecific location, and re-employed together with the delay blastinginitiation time series of the blast in order to predict time series dataof the ground vibration or noise of a single-hole blast of thesubsequent blast.

According to the blasting method of the present invention, the groundvibration or noise generated at a particular location upon a delayblasting can be controlled to a minimum without monitoring dominantfrequency of the ground and a waveform of a single-hole blast at alocation where ground vibration or noise becomes problematical prior toevery blast.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a waveform of vertical ground vibration at Location A. Thewave is produced by initiating two primers placed in water so as to haveinitiation timings of 10 ms and 40 ms, respectively, (i.e., aninitiation time interval of 30 ms). Each of the primers consists of anelectronic delay detonator and a water-gel explosive (100 g).

FIG. 2 shows a waveform of vertical ground vibration at Location A. Thewave is produced by initiating a primer placed in water so as to have aninitiation timing of 10 ms. The primer consists of an electronic delaydetonator and a water-gel explosive (100 g).

FIG. 3-1 shows a waveform of vertical ground vibration of a single-holeblast, which is predicted from the waveform shown in FIG. 1 according toa successive analytical prediction method described in the presentinvention.

FIG. 3-2 shows a waveform of vertical ground vibration of a single-holeblast, which is predicted from the waveform shown in FIG. 1 according tothe Fourier Transform described in the present invention.

FIG. 3-3 shows a waveform of vertical ground vibration of a single-holeblast, which is predicted from the waveform shown in FIG. 1 according tothe de-convolution method of the present invention.

FIG. 4-1 shows a computed waveform of vertical ground vibration atLocation A when, using the waveform of FIG. 3-1, a two-period delayblast with an initiation interval of 120 ms is effected according to thelinear superposition theorem.

FIG. 4-2 shows a computed waveform of vertical ground vibration atLocation A when, using the waveform of FIG. 3-2, a two-period delayblast with an initiation interval of 120 ms is effected according to thelinear superposition theorem.

FIG. 4-3 shows a computed waveform of vertical ground vibration atLocation A when, using the waveform of FIG. 3-3, a two-period delayblast with an initiation interval of 120 ms is effected according to thelinear superposition theorem.

FIG. 5 shows a waveform of vertical ground vibration at Location A. Thewave is produced by initiating two primers placed in water so as to haveinitiation timings of 10 ms and 130 ms, respectively (i.e., aninitiation interval of 120 ms). Each of the primers consists of anelectronic delay detonator and a water-gel explosive (100 g).

FIG. 6 shows a waveform of vertical grounds vibration at Location A. Thewave is produced by initiating five primers placed in water so as tohave initiation timings of 10 ms, 40 ms, 70 ms, 100 ms and 130 ms,respectively (i.e., initiation intervals of 30 ms). Each of the primersconsists of an electronic delay detonator and a water-gel explosive (100g).

FIG. 7-1 shows a waveform of vertical ground vibration of a single-holeblast, which is predicted from the waveform shown in FIG. 6 according toa successive analytical prediction method described in the presentinvention.

FIG. 7-2 shows a waveform of vertical ground vibration of a single-holeblast, which is predicted from the waveform shown in FIG. 6 according tothe Fourier Transform described in the present invention.

FIG. 7-3 shows a waveform of vertical ground vibration of a single-holeblast, which is predicted from the waveform shown in FIG. 6 according tothe de-convolution method of the present invention.

FIG. 8-1 shows a computed waveform of vertical ground vibration atLocation A when, using the waveform of FIG. 7-1, a five-period delayblast with an initiation interval of 90 ms is effected according to thelinear superposition theorem.

FIG. 8-2 shows a computed waveform of vertical ground vibration atLocation A when, using the waveform of FIG. 7-2, a five-period delayblast with an initiation interval of 90 ms is effected according to thelinear superposition theorem.

FIG. 8-3 shows a computed waveform of vertical ground vibration atLocation A when, using the waveform of FIG. 7-3 a five-period delayblast with an initiation interval of 90 ms is effected according to thelinear superposition theorem.

FIG. 9 shows a waveform of vertical ground vibration at Location A. Thewave is produce by initiating five primers which are placed in water soas to have initiation timings of 10 ms, 100 ms, 190 ms, 280 ms and 370ms, respectively (i.e., initiation intervals of 90 ms). Each of theprimers consists of an electronic delay detonator and a water-gelexplosive (100 g).

BEST MODE FOR CARRYING OUT THE INVENTION

Hereinafter, the blasting method of the present invention is illustratedin more detail with reference to Examples.

A plurality of primers, each of which consisted of an electronic delaydetonator (trade name: EDD) with an opportunity set initiation timingand a water-gel explosive (100 g) (trade name: Sunvex), was placed at adepth of 2 meters near the center of a pond (longer side: 25 m, shorterside: 25 m, depth: 4 m) so that the distance of each primer may be aboutone meter, and then initiated. The ground vibration (normal direction X,tangent direction Y, vertical direction Z) was monitored at a location100 meters away from the pond (hereinafter referred to as Location A) toconfirm the effects of the present invention.

EXAMPLE 1

Two electronic delay detonators, whose initiation timings were set so asto be 10 ms and 40 ms, respectively (i.e., an initiation interval of 30ms), were arranged individually in a water-gel explosive (100 g) andplaced in water. The detonators were exploded to monitor groundvibration thereby produced at Location A. Among the monitored waveforms,the one derived from the vertical ground vibration is shown in FIG. 1.An. electronic delay detonator, whose initiation timing was set so as tobe 10 ms, was arranged in a water-gel explosive (100 g) and placed inwater. The detonator was exploded to monitor ground vibration therebyproduced at Location A. Among the monitored waveforms, the one derivedfrom the vertical ground vibration is shown in FIG. 2.

From the waveform shown in FIG. 1, a vertical waveform of a single-holeblast producing the waveform of FIG. 1 was predicted. The waveformsobtained by the successive analytical prediction method, FourierTransform method and de-convolution method of the present invention areshown in FIGS. 3-1, 3-2 and 3-3, respectively.

Based on the linear superposition theorem, vertical waveforms of thesubsequent blasts of two-period delay blasts, whose initiation intervalswere set variously, were predicted using the above-predicted waveforms(FIGS. 3-1, 3-2 and 3-3). As a result, the maximum displacement velocityamplitude of the vertical wave at Location A was minimized when theinitiation interval was set at 120 ms. The predicted vertical waveformsof a two-period delay blast with an initiation interval of 120 ms, whichwere obtained according to the successive analytical prediction method,Fourier Transform method and de-convolution method of the presentinvention, are shown in FIGS. 4-1, 4-2 and 4-3, respectively.

In view of the above prediction, two electronic delay detonators, whoseinitiation timings were set at 10 ms and 130 ms, respectively (i.e., aninitiation interval of 120 ms), were arranged individually in awater-gel explosive (100 g) and placed in water. The detonators wereexploded to monitor ground vibration thereby produced at Location A.Among the monitored waveforms, the one derived from the vertical groundvibration is shown in FIG. 5.

Among the thus-obtained nine kinds of waveforms, the waveform shown inFIG. 2, which was derived from a single-hole blast, and the waveforms ofa single-hole blast shown in FIGS. 3-1, 3-2 and 3-3, which werepredicted according to the successive analytical prediction method, theFourier Transform method and the de-convolution method, were compared.As a result, it was found that these waveforms were very similar and thesuccessive analytical prediction method, the Fourier Transform methodand the de-convolution method were all advantageous in predictingwaveforms derived from a two-period delay blast. When the similarity ofthese waveforms was evaluated according to cross-correlationcoefficient, the correlation coefficients of FIG. 2 and FIGS. 3-1, 3-2and 3-3 were 0.88, 0.93 and 0.96, respectively. These results mean thatthe waveforms are similar in quantity, too.

Comparing the waveforms of a two-period delay blast shown in FIGS. 4-1,4-2 and 4-3, which were predicted at Location A based on the linearsuperposition theorem when a two-period delay blast was exploded with aninitiation interval of 120 ms using the waveforms of a single-hole blastpredicted according to the successive analytical prediction method, theFourier Transform method and the de-convolution method, with thewaveform of the vertical ground vibration shown in FIG. 5, thosewaveforms also very much resembled each other. The correlationcoefficients of FIGS. 4-1, 4-2 and 4-3 and FIG. 5 were 0.92, 0.92 and0.91, respectively.

EXAMPLE 2

Five electronic delay detonators, whose initiation timings were set soas to be 10 ms, 40 ms, 70 ms, 100 ms and 130 ms, respectively (i.e., aninitiation interval of 30 ms), were arranged individually in a water-gelexplosive (100 g) and placed in water. The detonators were exploded tomonitor ground vibration thereby produced at Location A. Among themonitored waveforms, the one derived from the vertical ground vibrationis shown in FIG. 6.

From the waveform shown in FIG. 6, a vertical waveform of a single-holeblast producing the waveform of FIG. 6 was predicted. The waveformsobtained by the successive analytical prediction method, FourierTransform method and de-convolution method of the present invention areshown in FIGS. 7-1, 7-2 and 7-3, respectively.

Based on the linear superposition theorem, vertical waveforms of thesubsequent blasts of five-period delay blasts, whose initiationintervals were set variously, were predicted using the above-predictedwaveforms (FIGS. 7-1, 7-2 and 7-3). As a result, the maximumdisplacement velocity amplitude of the vertical wave at Location A wasminimized when the initiation interval was set at 90 ms. The predictedvertical waveforms of a five-period delay blast with an initiationinterval of 90 ms, which were obtained according to the successiveanalytical prediction method, Fourier Transform method andde-convolution method of the present invention, are shown in FIGS. 8-1,8-2 and 8-3, respectively.

In view of the above prediction, five electronic delay detonators, whoseinitiation timings were set at 10 ms, 100 ms, 190 ms, 280 ms and 370 ms,respectively (i.e., an initiation interval of 90 ms), were arrangedindividually in a water-gel explosive (100 g) and placed in water. Thedetonators were exploded to monitor ground vibration thereby produced atLocation A. Among the monitored waveforms, the one derived from thevertical ground vibration is shown in FIG. 9.

The waveform shown in FIG. 2, which was derived from a single-holeblast, was compared with the waveforms shown in FIGS. 7-1, 7-2 and 7-3,which were predicted according to the successive analytical predictionmethod, the Fourier Transform method and the de-convolution method. As aresult, it was found that the waveforms very much resembled each otheras well as the comparison with those derived from a five-period delayblast. This means that the successive analytical prediction method, theFourier Transform method and the de-convolution method are always usefulto predict a waveform of a single-hole blast. The correlationcoefficients of FIGS. 7-1, 7-2 and 7-3 and FIG. 2 were 0.92, 0.96 and0.93, respectively.

Comparing the waveforms of a five-period delay blast shown in FIGS. 8-1,8-2 and 8-3, which were predicted at Location A based on the linearsuperposition theorem when a five-period delay blast was exploded withan initiation interval of 90 ms using the waveforms of a single-holeblast predicted according to the successive analytical predictionmethod, the Fourier Transform method and the de-convolution method, withthe waveform of the vertical ground vibration shown in FIG. 9, thosewaveforms also very much resembled each other. The correlationcoefficients of FIGS. 8-1, 8-2 and 8-3 and FIG. 9 were 0.86, 0.90 and0.89, respectively.

Industrial Applicability

The blasting method of the present invention is useful to reduce theground vibration and noise generated upon blasting.

What is claimed is:
 1. A blasting method which comprises conducting adelay blast at a particular location; predicting time series data of awaveform of ground vibration or noise at a remote location to begenerated by a hypothetical single-hole blast at the particular locationusing at least one of previous time series data of a waveform of groundvibration or noise generated by said delay blast and actually monitoredat the remote location, and the corresponding previous actually appliedinitiation time series of said delay blast; computing a delay blastinginitiation time series for a delay blasting, which provides a waveformof ground vibration or noise satisfying specific conditions, based onthe above-predicted time series data of a single-hole blast; andcarrying out a subsequent delay blast according to the computed delayblasting initiation time series.
 2. A blasting method according to claim1, wherein time series data of a waveform of ground vibration or noiseat a remote location to be generated by a hypothetical single-hole blastat the particular location is predicted by conducting a delay blast atthe particular location; then computing the Fourier Transform of thetime series data of a waveform of ground vibration or noise generated bysaid delay blast and actually monitored at the remote location, and thecorresponding actually applied initiation time series data of said delayblast to obtain corresponding spectrums; performing with the spectrums;and computing the Inverse Fourier Transform of the performed spectrum.3. A blasting method according to claim 1, wherein time series data of awaveform of ground vibration or noise at a remote location to begenerated by a hypothetical single-hole blast at the particular locationis predicted by conducting a delay blast at the particular location; andthen computing the cross-correlation sequence of time series data of awaveform of ground vibration or noise generated by said delay blast andactually monitored at a remote location, and the auto-correlationsequence of the corresponding actually applied initiation time seriesdata of said delay blast.